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1 .  INTRODUCTION 


In  the  investigation  of  the  vibrational  response  of  a  mortar  tube 
system  which  is  subjected  to  initial  conditions,  boundary  constraints 
and  applied  actions  with  the  normal  mode  method  of  dynamic  analysis  a 
nonlinear  equation  occurred  for  which  no  explicit  solution  could  be 
found  in  the  literature.  Since  this  equation  of  motion  describing  the 
vibrational  displacement  of  a  conical  ballistic  tube  with  a  concentric 
cylindric  bore  is  inherent  to  many  mortar  as  well  as  gun  tube  systems, 
an  explicit  analytical  solution  is  being  derived  for  a  variable  domain 
characteristic  to  ballistic  tubes. 

2.  FORMULATION  OF  THE  EQUATION  FOR  THE  TRANSVERSE  VIBRATIONS 
OF  A  CONICAL  BALLISTIC  TUBE  WITH  A  CYLINDRIC  BORE 

The  transverse  vibrational  displacement  of  a  beam  with  arbitrary 
cross  section  under  the  influence  of  applied  forces  is  governed  by  the 
equation  of  motion 


[El  y"(x,t)]"  +  pA  y(x,t)  =  F (x,t) , 


(2.1) 


where  y  is  the  lateral  displacement,  x  the  coordinate  measured  along  the 
axis  of  the  beam,  t  the  time  variable,  E  Young's  modulus  of  elasticity, 

I  the  moment  of  inertia  of  the  cross-sectional  area  with  respect  to  x, 
p  the  density,  A  the  cross-sectional  area,  and  F(x,t)  the  acting  force. 
The  symbols  '  and  •  denote  the  partial  derivatives 


,  respectively. 


Adopting  the  normal  mode  method  of  dynamic  analysis  the  solution  of 
the  inhomogeneous  equation  for  a  given  set  of  boundary  conditions  can  be 
approximated  to  any  degree  of  accuracy  by  superimposition  of  orthonormal 
functions  representing  eigen- solutions  of  the  homogeneous  equation  to 
the  specific  boundary  problems1.  Hence  the  problem  of  deriving  the 
solution  to  Eq.  (2.1)  for  a  set  of  boundary  conditions  is  reduced  to 
solving  the  homogeneous  partial  differential  equation  for  the  same  set 
of  boundary  conditions. 

[El  y"(x,t)]"  +  pA  y(x,t)  =  0  (2.2) 


*S.  Timoshenkoj  D.  H.  Young ,  W.  Weaver ,  Vibration  Problems  in  Engineering^ 
1974 ,  G.  Wiley  and  Sonsa  Ino.t  NY. 
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Since  flexural  rigidity  (El)  and  mass  unit  length  (pA)  are  time 
independent  for  practically  all  gun  tube  vibrational  problems,  we  can 
separate  the  variables  by  setting 


y(x,t)  =  g(x)  sin(at+a) , 


(2.3) 


where  a  is  the  frequency  of  vibration  and  a  is  the  phase  angle,  and 
obtain  a  homogeneous  fourth-order  ordinary  differential  equation  for  the 
displacement  amplitude  g(x) : 


[Elg  "  ] ’ '  -  a2pAg  =  0 


(2.4) 


The  fundamental  system  of  solutions 
linearly  independent  solutions,  g^, 

for  the  construction  of  the  general 


consists,  in  this  case,  of  four 
g2,  gg>  and  g^,  which  must  be  found 

solution  to  Eq.  (2.4): 


4 


g  = 


c  g 

n°n 


1 


(2.5) 


The  constants  c^ ,  c^,  and  c^  must  be  determined  to  within  an 

arbitrary  constant  in  each  particular  case  from  the  boundary  conditions 
at  the  end  of  the  beam.  In  general,  the  end  boundary  conditions  for  gun 
tube  vibrational  response  problems  are  physical  constraints  on  deflection 
y,  slope  y',  bending  moment  [Ely  *  * ] ,  and  shearing  force  [Ely'']'.  For 
example,  at  a  fixed  end  the  deflection  and  slope  are  equal  to  zero 
(y  =  0,  y'  =0);  at  a  simply  supported  end  the  deflection  and  bending 
moment  are  equal  to  zero  (y '  =  0,  [El  y"]  =  0);  and  at  a  free  end  the 
bending  moment  and  the  shearing  force  both  vanish  ([Ely'']  =  0,  [Ely ' * ] 1 
=  0) .  At  an  interior  boundary  the  continuity  of  deflection,  slope, 
bending  moment  and  shearing  force  must  be  preserved.  Since  these  functions 
appear  in  the  formulation  of  the  eigenvalue  problem,  we  will  derive 
explicit  analytic  expressions  for  them  concurrently. 

In  this  study  we  are  only  interested  in  the  case  where  Young's 
modulus  of  elasticity  and  the  density  are  constant  and  the  beam  is  a 
conical  ballistic  tube  with  a  cylindric  bore.  The  geometry  of  the  beam 
is  outlined  in  Figure  2.1. 
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UP 


Substituting  Eqs.  (2.6),  (2.7),  and  (2.8)  into  Eq.  (2.4),  we  obtain 

&4’1)8"]  "  e7&2  (*2  ’ x)8  =  °*  (2-9) 

Rr  (Rl  Rr\  x  +  Lr 
-r  *\r  'r  I  \  +  Lr  * 
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3.  DERIVATION  OF  THE  FOUR  LINEAR  INDEPENDENT  SOLUTIONS 


The  explicit  form  of  Eq.  (2.9)  suggests  the  following  coordinate 
transformation: 


where  the  new  variable  z  is  now  dimensionless  and  defined  in  the  domain 

.  For  most  current  weapon  systems,  the  variable  z, 

which  is  the  ratio  of  wall  thickness  to  bore  radius  of  the  ballistic 
tube,  lies  within  the  doman  0  <z<  1.  Therefore,  the  investigation  will 
mainly  be  restricted  to  this  geometric  domain  of  interest.  Carrying  out 
the  mapping  of  Eq.  (2.9)  from  the  x-space  onto  the  z-space,  we  get 


uT(2  +  z)zg  =  0, 


(3.2) 


where  the  separation  constant  a  is  redefined  as 


a  = 


kw 


with  k 


1 


-  L*  + 


L  • 

r 


(3.3) 


The  ordinary  differential  Eq.  (3.2)  has  regular  singular  points  at  z  =  0, 
-2,  (-1  +  i),  (-1  -  i),  and  an  essential  singularity  at  z  =  «.  Since 
z  =  0  is  only  a  regular  singularity,  we  can  obtain  a  solution  in  the 
neighborhood  of  the  origin  by  a  series  expansion.  We  try 

oo 

g(z)  =  Ear,  ao  f  0,  (3.4) 

v=0 


with  the  exponent  n  and  all  the  coefficients  a  still  undetermined.  By 
differentiating  twice,  we  obtain 

co 

8zz  =  av(n  +  vHn  +  V  -  l)zn+v"2.  (3.5) 

v=0 
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By  multiplying  this  expression  by  zm  and  differentiating  it  again  twice, 

we  get 

00 

pSzzjzz  =  av  (n  +  v)(n  +  v  •  1)(n  +  m  +  v  "  2)  (3.6) 

^  J  v=0 


x(n  +  m  +  v  -  3)z(n+m+v-4^. 


By  substituting  into  Eq.  (3.2),  we  have 

oo 

'y  ^  (n  +  v)  (n  +  v  -  1)  (n  +  v  +  2)  (n  +  v  +  l)zn+v 
v=0 

+4(n  +  v  +  1) (n  +  v)zn+v'1  +  6(n  +  v) (n  +  v  -  l)zn+v ~2 

+4(n  +  v  -  l)(n  +  v  -  2)zn+V_3]  -  u>4(2  +  z)z  ax)zn+V  =  0.  (3.7) 


The  uniqueness  theorem  of  power  series  requires  that  the  coefficients  of 
each  power  of  z  on  the  left  hand  of  Eq.  (3.7)  must  vanish  individually. 

The  lowest  power  of  z  appearing  in  Eq.  (3.7)  is  zn”3  for  v  =  0. 

The  requirement  that  the  coefficient  is  zero  yields 

4n(n  -  l)*(n  -  2)aQ  =  0.  (3.8) 


Since  aQ,  by  definition,  is  the  coefficient  of  the  lowest  non-vanishing 
terms  of  the  series,  we  have 


n(n  -  l)2 (n  -  2)  =  0. 


(3.8.1) 


The  roots  of  this  indicial  equation  are  n  =  0,  1,  1,  2  resulting  in 
three  independent  power  series 

OO  00  00 

gl  ■  2/lv2'’  ’  h  "  1  S  a2vZ”  •  *3  s  z2  2Z  a3vz'’-  f3-9> 

v=0  V=0  v=0 
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Since  n  =  1  is  a  double  root,  a  fourth  linearly  independent  solution  may 
be  obtained  by  multiplying  the  regular  solution  g-j(z)  by  a  logarithmic 
term  and  adding  another  power  series: 

00 

g4(z)  =  g2Cz)ln  2  +  2a4vzV+nr  (3.10) 

v=0 


For  the  eigenvalue  problem  we  need,  in  addition  to  the  deflection  g(z), 
the  slope  g^,  the  bending  moment  [EIgzz]>  and  the  shearing  force  [EIgZ7_l7- 

To  directly  obtain  explicit  expressions  for  these  physical  functions,  we 
will  slightly  modify  the  approach.  Integration  of  Eq.  (3.2)  leads  to 


|(z3  +  4z2  +  6z  +  4^zgzzJz  ®  Kj  +  a)  •/  dzz (2  +  z)g(z) 

_  r  2  z  1  kn+2+v 

^^anv  Ik  +  2  +  v  +  k  +  3  +  v  I  Z 

o  L  n  n  J 


„n  4 

+  (JL> 


(3.11) 


(3.11.1) 


...for  n  =  1,  2,  3  and  k  =  0,  1,  2. 
'  n  ’ 


„4  4  ,  , 

=  Kj  +  a)  \  In  z 


{  1  [*2  +^2  k2  -] z  2 


2v  |  (  k?  +  2  +  V 


■7  (k2  *  3  *  ") 


z  1  k2* 

7777FJZ 


2+v 


4v 


_ 2 _  z  1  k4+2+V  | 

+  2  +  v  +kq  +  3  +  vJ  / 


(3.11.2) 


for  n  =  4  and  k,  to  be  determined. 

-4 


Integrating  the  above  equation  again,  we  have 


( 


3  2 

z  +  4z  +  6z  +  4 


')z«zZ  *  Ko  *  K1  1  ‘  "V  dZ  J 


dz’z'  (2  +  z')  g(z’) 
(3.12) 
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Substituting  the  power  series  Eqs.  (3.9)  and  (3.10)  into  this  equation,  we 
have 


£anv(kn  *  v)  (kn  *  v  -  *)  (! 


3  2 

z  +  4z  +  6z  +  4 


„n  „n  4  ^ 

K  +  K.  z  +  a)  >  a 
o  1  /  j  nv  k 


00 

v,  i  r  2 

Z^“nv  k  +  3  +  v  k  +2  +  v 
n  n  L  n 


1  k 

z  n 

- a -  z 

<  +  4  +  v 

n 


0 

k  +v+3 


(3.12.1) 


for  n  =  1,  2,  3  and  kn  =  0,  1,  2; 


In  z 


Sa2v(k2  +  V)  (k2  +  V  '  X)  (z3  +  4z2  +  6z  +  4)z  2 

n  '  ' 


Z{aav  [2(k2  **)  -  "4  *  v(k4  ‘  “)  (k4  -  -  >)} 


,  ~  \  k.+v-l 

x  |z~  +  4z"  +  6z  +  4lz 


(2 


ko  *  K1  2  *  I1"  2  2_.a 2v  k 


00 

I].a2v  k2  ♦  3  .  v  [k2  *  2  .  v 

z  1  k2+V'3  vM  1  f _ 2 

k^  +  4  +  vjZ  +  ^  a2v  k2  +  3  +  v  |  k2  +  2  + 

(k2  ♦  2  *  v  *  k2  *  3  +  v  )  *  k2  *  4  ♦  v  ^k2  *  3  -  v 


+  k2  +  4  +  v)j 

k4  +  4  +  j 


k  -k 

2  K4  1 

z  +  a 


4v  k,  +  3  +  v 
4 


[k4  * 


k  .+v+3 
4 


(3.12.2) 
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for  n  =  4  and  k„  to  be  determined. 

4 

The  uniqueness  of  power  series  requires  that  the  total  coefficient  of 
each  power  of  z  vanishes  all  by  itself.  For  the  indicial  equation  roots 
n  =  1 ,  2 ,  and  3 ,  we  have 

z°:  -kJ  =  0;  =  0:  -Kq  =  o 


Kq  =  0;  n  =  1,  2,  3 


1 

z  : 


-kJ  +  2.1.4a12  =  0;  -kJ  +  2.1.4a21  =  0; 


■K1  +  2 ' 1 *4a30  "  0 


K1  8an,3-n;  n  "  1,2,3 

2 

z  :  2.1.6a^2  +  3.2.4a^3  =  0 

2.1.6a2^  +  3.2.4a22  =  0 
2.1.6a3Q  +  3.2.4a31  =  0 


a  a  =  -  a  -  ;  n=i,2,3 

n,4-n  2  n^-n*  *  * 


3 

z  : 


03 


2.1.4a12  +  3.2.6a13  +  4.3.4aH  3>2  .  —10 


.  2a. „  =  0 


2.1.4a2^  +  3.2.6a22  +  4.3.4a23  =  0 


2.1.4a3Q  +  3.2.6a3^  +  4.3.4a32  =  0 


(hquations  continued  on  next  page) 
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,Za13  +  /,ial2  "  3.2  U 


4 
z  : 


i  r  3 

*  14  4.3  [2  * 

an,5-n  =  '  471  |j  *  3,2an,4-n  +  2*lan,3-nJ  ;  n=2’3 

"  473  (2all  +  alo)  = 


2*lal2  +  3.2.4a^3  +  4.3.6a^4  +  5.4.4a^ 


w 


2.1a21  +  3.2.4a22  +  4.3.6a23  +  5;4.4a24  -  ~j  .  2&2Q  =  0 


2. la^Q  ^  3.2.4a3^  +  4.3.6a32  +  5.4.4a33  =  0 


'15 


3aH  +  3.2a13  +  4  .  2.1a12  - 


24 


3a23  +  3.2a22  +  4  .  2.1a2l  -  -~ 


33 


.  J_  ("l  4 

5.4  |  2  ‘ 

•  0K-.  •  -,)] 

-  JL  |l  4 

5.4  [2  *  4- 

« (r)  2-2o] 

-  573-  |y  .  4.3a32  ♦  3.2a31  ♦  j  .  Z.laJ 

74  (2a12  *  all)  ’ 

"  574  (2a21  *  a2o)  =  0 


3.2a13  +  4.3.4a14  +  5.4.6a15  +  6.5.4a16  -  \2&u  +  an  /  =  0 


3. 2^22  ^  ,^,^a23  6.5.43.2^ 


3.2a31  +  4.3.4a32  ♦  5.4.6a33  +  6.5.4a34  -  ^  .  2a3Q  =  0 


^  n ,  7  -n  6 
(Equations  continued  on  next  page) 


1  I3  r  , 

F"  It".  5 . 4a  , 

.5  |2  n,6-n 


+  4.3a  r  +  j-  .  3.2a 

n,5-n  4  n,4-n 
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34 


5.4  ^4  )  (^an,3-n  +  an,2-n^ 

— r*  s 

6.5  [2  b' 

(£)  2a3o] 


;  n  =  1,  2 


4a33  +  4,3a32  +  4  •  3,2a31  "  5^4 


6 

z  : 


4.3a14  +  5.4.4alc.  +  6.5.6alft  +  7.6.4a17 


n,8-n 


-  (2ai3 +  ai2) =  0 

'  6  75  ^2a22  +  a21^  =  0 

‘  675  (2a31  +  a3o)  =  0 

Y  .  4 . 3a  c  -  ( 7-  J  ( 2a  „  +  a  ,  ^ 

4  n,5-n  6.5  \4  /  \  n,4-n  n,3-nJ 


4.3a?1.  +  5.4.4a,,4  +  6.5.6a?,-  +  7.6.4a?A 


4.3a__  +  5.4.4a__  +  6.5.6a_.  +  7.6.4a_r 
32  33  34  35 


0  .  6.5a  _  +  5.4a  , 

2  n,7-n  n,6-n 


n  =  1,  2,  3 


.  .  .  (3.13) 

Without  loss  of  generality  we  can  set  the  coefficients  a^,  a^»  a^,  a21> 
&22’  anc*  a23  t0  zer0,  By  inspection  and  mathematical  induction  we  obtain 
the  following  recurrence  relations  for  the  coefficients: 


(Equations  continued  on  next  page) 
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a20  *  °’ 

a21 

CM 

CM 

n 

a30  *  °> 

a31 

1 

a  -  — 

2 

1 

an,v+l 

(n  + 

v)  (n 

24 


aT0  =  TA  *  a 


120  \4  / 


32  24  *  33  16 


34 


=  _L_  +  _i_  41) 

240  300  \4  / 


(n  +  v  -  1)  (n  +  v  -  2) a 


n,v 


+  (n  +  v  -  2)(n  +  v  -  3)an  ^  (n  +  v  -  3)  (n  +  v  -  4) 


x  a 


n,v-2  "  (n  +  v  -  1)  (n  +  v  -  2) 
..  for  n  =  1,  2,  3  and  v  =  4,  5,  6,  .... 


^4  )  ^2an,v-3  +  an,v-4jJ 


(3.14) 


To  derive  the  irregular  solution,  we  apply  the  uniqueness  theorem  of  power 
series  to  Eq.  (3.12.2).  Because  the  logarithmic  terms  vanish,  we  have 

z°:  "4 * la20  *  K0  =  0 


K0  4a20 


z  :  -6.1aon  -  4.3a.,,  +  K,  =  0 


K1  =  6a20 


z  :  -4.1a2Q  -  6.3a2l  -  4.5a22 


This  expression  is  ^  0;  to  cancel  it  we  must  require  that 


k.  =  3. 
4 


-4,la20  "  6,3a2i  "  4,5a22  "  4,3,2a40  "  0 


(Equations  continued  on  next  page) 
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a40  "6  a20 

-l.la2o  “  4,3a21  "  6,5a22  "  4,7a23  "  6'3'2a40  "  4,4,3a4i  =  0 
a41  =  48  a20 


4  7  us 

z  :  -3a21  -  4.5a22  -  6.7a23  -  4.9a24  -  2  2a20  '  4,3,2a40 

.  4  J 


5 

z  : 


-6.4. 3a, j  -  4.5.4a42  =  0 


42 


[— *  —  6^)1 
| 160  3600  \4  /  | 
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-5a22  -  4.7a23  -  6.9a24  -  4.11a25 


9w 


(4.5) 


2  (2a21  +  a2o) 


6 
z  : 


-3.2a40  -  4.4.3a41  -  6.5.4a42  -  4.6.5a43  =  0 


43  6 


.5  ’  (  |_2  *  5 ' 4a42  +  4 ' 3a4 1  +  4  ‘  3,2a4oJ 

lla25  +  2  ‘  9a24  +  7a23  +  4  *  5a22  +  ~~I 
L  (4.5 

x  (r)  (2a2i  *  a2o)]  ) 


7a23  "  4 ' 9a24  '  6,lla25  '  4 * 13a26 


lid) 


(5.6) 


2  {^a22  +  a21^ 


-4.3a41  -  4.5. 4a, 2  -  6.6.5a43  -  4.7.6a,,  ♦  ^  .  2a, Q  =  0 


(Equations  continued  on  next  page) 
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44 


7'6  i[2 

-  67s  (f)  2a4„]  *  [ 


6.5a.  _  +  5.4a...  +  -r  .  4.3a., 
43  42  4  41 


,  1  11 

*7'7a»*^7 


13a26  +  2  •  lla25  +  9a24 

(?“)  (2a22  +  a2l)Jj 


"9a24  "  4,lla25  "  6,13a26  "  4,15a27 


13  0) 


(6.7) 


2  (2a23  +  a22) 


-S.4a42  -  4.6.5a43  -  6.7.6a44  -  4.8.7a45  *  S—  (; 


2a41  *  a40 


')  ■“ 


45 


.  7 

8.7  ||2  • 


^a44  +  ^"^a43  +  4  *  3,4a42 


2a41  +  a40 


+  lla25  +  4  •  9a24  +  2 

4  4  (6. 7j 


,)|  *  |15a27  *  I 


.  13a 
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2a23  +  a22 , 


(3.15) 


From  this  pattern  we  can  deduce  the  general  recurrence  formula: 


4 , v+1  '  (4 


1  [U 

*  v)(3  ♦  v)  U2 


*  a4 ,v-l  *  ?(1  *  v>v  a4,v-2  -  U 


j  (3  +  v)(2  +  v)a4  v  +  (2  +  v) (1  +  v) 

_ l _ 

+  v)(2  +  v)  Y  J 
y(2a4,v-3  *  a4,v-4)|+  *  7>  a2,v+3  *  f  <2”  *  5’a2,v.2 

* (2u  +  3)  a2-  *  *  ^  ^ 2) 2  ¥) 

*  K"-l  *  a2.a-2)]  j 


(3.16) 
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To  ascertain  that  Frobenius'  method  yields  series  solutions  which  not 
only  satisfy  our  fourth-order  differential  equation  but  also  converge 
over  the  region  of  interest  (0  <z<_  1) ,  we  have  to  determine  the  convergence 
of  the  series.  The  recurrence  formulae  (3.14)  and  (3.16)  exhibit  a 

4\ 

dependence  for  the  coefficients  a  which  is  of  the  following 


structure : 


an0’  anl * 
an4 ’  an5’ 
an8’  an9’ 


(3.17) 


Therefore,  we  will  apply  the  D'Alembert  ratio  test  to  the  ratio  of  the 
coefficients  an  and  a^  ^ .  For  large  v,  we  can  approximate  the 

recurrence  formulae  for  the  regular  solutions,  Eq.  (3.14),  and  for  the 
nonlogarithmic  part  of  the  irregular  solution,  Eq.  (3.5),  by 


a 


n,  v+4 


3 

2  an,v+3 


an,v+2 


i_ 

4  an,v+l 


(3.18) 


In  order  to  deduce  a  recurrence  formula  which  contains  only  terms  propor¬ 
tional  to  a  , ,  a  ,  a  .,...we  will  use  the  recurrence  formula  for 
n,v+4  n,v  n,v-4 

the  (v+3)-th  term  in  addition  to  that  of  (v+4)-th  term. 


n,  v+3 


3 

2  an,v+2 


a  , 
n,  v+1 


1 

~r  a  + 
4  n,  v 


(3.19) 


Addition  and  subtraction  of  Eqs.  (3.17)  and  (3.18)  yield 

an,v+3  5  an,v+4  2  an,v+2  +  20  an,v  +  ^  (v)  (3.20) 

dn,v+l  5  an,v+4  an,v+2  10  an,v  +  \y)  '  (3.21) 
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These  relations  still  contain  the  coefficient  a  -.To  eliminate  the 

n,v+2 

an  v+2  dePendency>  we  aPPly  relations  (3.20)  and  (3.21)  to  the 
[ (v+3) -2, (v+1) -2]  and  [(v+3)-4,(v+l)-4]-th  coefficients. 


a  , 
n,v+l 


4 

5  an,v+2 


+  20  an, v-2 


a 

n,v-i 


4 

5  an,v+2 


a 


n,v 


3_ 

10  an,v-2 


a 


n,  v-1 


4 

a 

b  n,v 


1_ 

2  an,v-2 


+ 


1_ 

20  an,v-4 


an,v-3 


4 

-=-  a 

5  n,v 


an,v-2 


(3.22) 

(3.23) 

(3.24) 

(3.25) 


By  successive  substitution  of  Eqs.  (3.21),  (3.22),  (3.23),  and  (3.24) 
and  simple  manipulation  we  obtain  the  following  identity  in  the  form  of 
a  limit: 


, -  f  an,v+4 

limJ  —  • 

v-*»  (  n,v 


an,v  +  _3  an,  v 

a  .  16  a  , 

n,v-4  n,v-4 


(3.26) 


Setting 


R  =  lim 


( 


a  a 
n,v-*-4 

a 

n,v 


lim 


( 


an,v-4  / 


(3.27) 


and  taking  the  limit  as  v-**>,  we  derive  a  quadratic  equation  for  the 
ratio  of  the  coefficients: 


R 


2 


1 

64 


0  . 


(3.28) 


Its  roots  are 


D  1  ,  1 

R  =  -  —  and  — 
4  lb 


(3.29) 
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Using  these  values,  we  can  establish  a  lower  bound  of  the  radius  of 
convergence  for  the  series  solution  as 


r*  =  min 


1_ 

4  ’ 


min  (/2,2)  = 


(3.30) 


which  is  equal  to  the  distance  between  the  point  z  =  0  and  the  regular 
singular  points  z  =  (-  1  ±  i) .  Since  our  power  series,  Eqs  (3.9)  and 
(3.10)  are  of  the  form 


gn(z) 


n- 


T 


n ,  v 


+  6 


4n 


oo 


(z  In  z) 


(3.31) 


(n  =  1,  2,  3,  4;  6^ . .  .Kronecker  delta  function); 


they  converge  at  the  point  z  =  0.  The  power  series  solution  to  Eq.  (3.2), 


g(z) 


4 

2Yn8nt2)-  CV" 
1 


arbitrary  constants) , 


(3.32) 


-j  Q 

is  convergent  for  z  =  re",  0  £  r  <  r*,  -  tt  <  0  <_  tt,  in  the  z-plane  cut 
along  the  negative  real  axis  and  uniformly  and  absolutely  convergent  for 
any  interior  domain,  r  <_  r**,  where  0  <_  r**  <  r*  =  /2  . 

•  y 


24 


Since  each  of  the  terms  g  (z)  =  a  zn+v  +6.  (z  In  z)  a„  zn, 

nv  nv  4n  2v 

n=  1,  2,  3,  4,  is  a  continuous  function  of  z  and  the  partial  solutions 
gn>  n  =  1,  2,  3,  4,  converge  uniformly  for  |z|  the  power  series 

solution  g(z),  Eq.  (3.31),  is  a  continuous  function  in  our  domain  of 
uniform  convergence.  With  g^U)  continuous  and  g(z)  uniformly  convergent, 

both  the  differentiated  and  the  integrated  series  are  power  series  which 
again  are  continuous  functions  and  have  the  same  radius  of  convergence 
as  the  original  series. 

By  use  of  Eqs.  (3.31),  (3.16),  and  (3.14)  and  performing  some  trivial 
mathematical  manipulations  we  can  express  deflection,  slope,  bending 
moment  and  shearing  force  as 

4 

gO;  z)  =  2^jyk k;  1  =  l>  2>  3*  4»  (3.33) 

k=l 


where  the  parameter  £  refers  to  the  following  physical  functions: 


£  =  1  ...  g(z);  deflection 
£  =  2  ...  g(z) J  slope 


£  =  3  . 


£  =  4  . 


((z  +  l)4  -  1^  — j  g(z)J  ;  modified  bending  moment 
4t-  l  (z  +  l)4  -  li  =-g(z)i  ;  modified  shearing  force. 

LV  7  dz2  J  (3. 


34) 


The  Yj,  are  the  coefficients  of  the  four  linearly  independent  solutions 

g(l»  k;  z),  k  =  1,  2,  3,  4,  and  must  be  determined  to  within  an  arbitrary 
constant  in  each  particular  case  from  the  boundary  conditions.  The  four 
linearly  independent  solutions  and  their  modified  derivatives  are  given 
by 
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g(£,  k;  z)  =  c(£,  k,  n)zn  +  <$4kg(*,  2;  z)ln  z,  (3.35) 


n=0 


where  the 

coefficients  c(£,  k. 

n) 

are  defined  as 

o 
/ — \ 

k, 

-n) 

=  0  for  all  n 

=  ( 

),  1,  2,  3 

9  •  •  • 

(3.36) 

c(l. 

1, 

0)  = 

1;  c(l,  2, 

0) 

=  c(l,  3, 

0)  = 

c(l. 

** 

O 

N - * 

II 

0 

c(l. 

1, 

1)  = 

0;  c(l,  2, 

1) 

=  l;  c(l,  3, 

1)  = 

c(l,  4, 

1) 

=  0 

c(l, 

1, 

2)  = 

c(l,  2,  2)  = 

0; 

c(l,  3, 

2)  = 

i; 

c(l,  4, 

2) 

=  0 

c(l, 

k. 

n)  = 

.  i  | 

n(n  -  1)  j 

[f 

(n  -  1)  (n 

-  2)c(l, 

k,  n  -  1) 

+ 

(n  -  2) 

*  (n  -  3) c (1 ,  k,  n  -  2)  +  J-  (n  -  3)  (n  -  4)c(l,  k,  n  -  3) 

(n  -  2)  \jr)  (2c(1>  k>  n  '  4)  +  c(i,  k,  n  -  5))| 


(n  -  1) (n 


[■ 


+  64k  | (2n  -  1)  c(l,  2,  n  -  1)  +  y  (2n  -  3) 


X  c(l,  2,  n  -  2)  +  (2n  -  5)  c(l,  2,  n  -  3)  +  j  (2n  -  7) 


x  c(l,  2,  n  -  4)  + 


2n  -  3 


(n  -  1) 2  (n  -  2) 2 


(2c (1,  k,  n  - 


5) 


c(l,  k,  n  -  6))]  j  ;  n  =  3,  4,  5,  ... 


(3.36.1) 
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c(2,  k,  n)  =  nc(l,  k,  n  +  1)  +  64kc(l,  2,  n  +  1);  for  all  n=0,l,2,3,. 

(3.36.2) 

c(3,  1,  0)  =  c(3,  2,  0)  =  c(3,  .3,  0)  =  0;  c(3,  4,  0)  =  4 
c (3,  1,  1)  =  c(3,  2,  2)  =  0;  c(3,  3,  2)  =  8;  c(3,  4,  2)  =  6 

c(3,  k,  n)  =  T)  |  [2c (1 ,  k,  n  -  3)  +  c(l,  k,  n  -  4)] 


"64k  Krn”-'~1)  [2c(1'  2»  n  -  3)  +  cd.  2,  n  -  4)]  | 


n  =  2,  3.  4.  ... 


f 3. 36. 31 


c (4 ,  1,  0)  =  c(4, 

/  4  \ 

2,  0)  =  0; 

( 

c(4,  3,  0)  =  8; 

c(4,  4 

r(A  1  n'l  -  f  1 

«  71  1 

w  >  1X  >  »*;  -  i  n  / 

|  9  ^ 9 

11  -  ^  J  T  f  IN.  , 

11  "  'J)\ 

h  [2c (1 ,  2, 

n  -  2)  +  c(l,  2, 

n  -  3)] 

n  = 

n 

12  3 

A  )  )  •  •  • 

(3.36.4) 


ana  o  .  is  tne  KronecKer  aelta  function  16  =0  for  m  f  n  and  6  =1 

*+k  \  \  mn  mn 

for  m  =  n  j . 

4.  FREE  LATERAL  VIBRATIONS  OF  A  CONICAL  BALLISTIC 
TUBE  WITH  A  CYLINDRIC  BORE 

The  determination  of  the  eigenvalues  of  the  separation  constant 
2 

a  =  kw  and  the  relative  values  of  the  coefficients  y^,  Y2>  Yj>  and  Y^ 

are  contingent  on  the  boundary  conditions  at  the  ends  of  the  beam.  Let 

4 

us  assume  that  the  solution  of  Eq.  (3.2),  g(l,  z)  =  ^Jy^  g|i,  k;  w4;z]  , 


is  constrained 
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(4.1.1) 


at  z  =  zr  by  gC^,  zy)  =  g(£2,  zr)  =  0,  ^  /  Z2>  and 
at  z  =  by  g (&3 ,  z£)  =  g(«-4 ,  z£)  =  0,  £3  t  (4.1.2) 


where  ^3*  and  ^4  are  any  °f  the  physical  functions  given  in  Eq. 

(3.34j.  In  matrix  formulation,  the  above  boundary  value  equations  can 
be  written  as 


8  y  1  \ 

8  (£2,1; 


l; 


,1; 


g 


g 


r  ’  w4) 

zr;  J ) 


V  “ 


V  w 


(4.2) 


A  non-trivial  solution  of  Eq.  (4.2)  exists  only  if  the  determinant  of 
the  coefficients  is  equal  to  zero 

g  ^  1  >  1  >  ^  ^  ^  j  g^^»3>  z^,  CO  J  g^tj,4,  z^,  a)  J 

i 

^2  5  4  *  V 

|s  [  ~  z *  *  >  z  0  y  ~  )  8  ^  >  -  5  -  q  >  ~  )  -  0  >  —  /  S  ^  3  i  ^  ~  ^ 


g!£9,l,  z  *  a)  )  g(l0$2;  z  ,  u)  I  gl£«,3,  z  ,  to  I  g 


g^4,!;  zA;  w4)  g^4,2;  z^;  u4)  g  ^4>3;  z^;  w4)  g^4,4;  z^;  w4^ 


=  0 


(4.3) 
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This  equation  furnishes  the  means  for  the  determination  of  the  eigenvalues 
(u4)  and,  consequently,  from  Eq.  (4.2)  the  ratios  of  the  coefficients 
Y2/Yi»  Y-j/Yj  and  y  /y  and  from  Eq.  (3.33)  the  eigenfunctions  gn(l;z)  as 
well . 


The  roots  of  the  characteristic  Eq.  (4.3)  represent  the  eigenvalues 
of  the  system  and  correspond  to  the  natural  frequencies  of  the  beam. 

Using  the  boundary  conditions,  Eq.  (4.1),  and  the  orthogonality  relation 

for  the  eigenfunctions,  one  can  show  that  the  eigenvalues  (“4)n  must  be 

positive  or  zero.^  Due  to  the  complexity  of  the  functions  g(£,  k;  z;  to4) 

the  eigenvalues  (“>4)n  cannot  be  determined  analytically  but  must  be  found 
numerically. 

As  a  numerical  example,  we  will  consider  the  eigenvalue  problem  of 
the  free  transverse  vibration  of  a  circular  conical  tube  with  a  circular 
cylindric  bore  where  the  bottom  of  the  cone  is  a  fixed  end  and  the  top 
is  allowed  to  vibrate  freely  (Figure  4.1)  and  compare  the  solution  to 
this  problem  with  the  free  transverse  vibrations  of  a  circular  tube. 


Figure  4.1.  Geometry  of  Numerical  Example 

.....  _  3.3 

We  select  steel  with  a  density  of  7.84  x  10  kg/m  and  a  Young's 

modulus  of  elasticity  of  2.1  x  loll  Pa  as  the  tube  material,  a  beam 

length  of  1.00m,  an  interior  radius  of  0.05m  and  for  the  circular  pipe 

an  exterior  radius  of  0.06m  and  determine  the  exterior  radii  R  and  R 

by  keeping  the  volume  constant  r 


'A.  S.  alder,  "tree  and  toraed  Vibrations  of  a  Tapered  Cantilever  Beam," 
University  of  Delaware;  M.  A.  Thesis,  June  1956. 
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Tables  4.1,  4.2,  and  4.3  are  representative  samples  of  the  obtained 
numerical  results. 


TABLE  4.1.  NATURAL  FREQUENCIES  OF  A  CONCENTRIC  CIRCULAR 
CONE  WITH  A  CYLINDRIC  BORE  AS  FUNCTION  OF  THE 


MODE  OF 

VIBRATION  FOR  R  = 

0.95 

n 

n 

(JL) 

A 

n 

n 

n 

i 

9.23175 

6.32619 

.45929 

2 

18.66643 

15.83683 

.17867 

3 

28.86741 

26.50023 

.08933 

4 

39.18171 

37.09654 

.05621 

5 

49.56180 

47.69570 

.03913 

6 

59.98144 

58.29474 

.02893 

7 

70.42839 

68.89379 

.02227 

8 

79.49283 
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TABLE  4.2.  NATURAL  FREQUENCIES  OF  A  CONCENTRIC  CIRCULAR 
CONE  WITH  A  CYLINDRIC  BORE  AS  FUNCTION  OF  TOE 
MODE  OF  VIBRATION  FOR  R  =  -0.95 


n 

0 

n 

0) 

n 

i 

3.67169 

6.32619 

-.41960 

2 

13.36686 

15.83683 

-.15596 

3 

24.81891 

26.50023 

-.06345 

4 

35.66632 

37.09654 

-.03855 

5 

46.41514 

47.69570 

-.02685 

6 

57.11542 

58.29474 

-.02023 

7 

68.89379 

8 

79.49283 

TABLE  4. 
CONE 

3.  NATURAL  FREQUENCIES  OF 
WITH  A  CYLINDRIC  BORE  AS 
RADII 

A  CONCENTRIC  CIRCULAR 
FUNCTION  OF  VARIOUS 

lRl 

AjC+lRl) 

“l 

fljC-iRl) 

6.01764 

90233 

6.66899 

,80417 

10.15110 

7.48786 

5.16764 

,70550 

8.52967 

60631 

12.37331 

9.92501 

7.75756 

,50661 

11.87835 

,40637 

17.10772 

14.80818 

12.69508 

,30560 

19.69097 

,20429 

31.65520 

29.45615 

27.34957 

,10243 

58.75086 
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Tables  4.1  and  4.2  contain  the  first  eight  natural  frequencies, 

of  a  concentric  circular  cone  with  a  concentric  circular  cylindric  bore 
for  two  selected  values  of  and  Rr  expressed  by  the  ratio  R  = 

^  -  R^J/max  ^R^  -  R^.  The  wn's  are  the  corresponding  eigenvalues  for 

the  circular  tube.  The  quantity  | /^expresses  the  deviation 

of  from  to  .  Table  4.3  exhibits  the  dependence  of  the  first  natural 

frequency  on  the  radii  of  the  cone. 

The  implication  of  the  numerical  result  to  the  geometrical  descript¬ 
ion  of  gun  tubes  for  lateral  tube  motion  can  best  be  seen  from  Figure 
4.2.  Here  we  plotted  the  deviation  of  the  natural  frequencies  of  the 
concentric  circular  cone  with  a  concentric  circular  cylindric  bore  from 
the  circular  tube  versus  the  normalized  difference  between  the  left  and 
right  radii  of  the  cone.  The  plot  displays  clearly  that  for  eigen- 
frequencies  in  the  lower  principal  mode  of  vibration  the  commonly  used 
approximation  of  a  conical  gun  tube  by  a  cylindric  beam  is  not  justified. 
However,  for  the  higher  principal  modes  of  lateral  vibration  we  can  well 
substitute  the  eigenvalues  of  the  beam  for  the  correct  ones. 

Having  the  natural  frequencies  determined  we  can  use  Eq.  (3.33)  to 
compute  the  deflections  of  the  centerline  of  the  concentric  circular 
cone  with  a  concentric  cylindric  bore  in  the  r-th  mode.  Figure  4.3  is  a 
representative  example.  The  values  on  the  abscissa  correspond  to  the 
deflections  of  the  beam  axis  normalized  by  the  maximum  displacement  and 
on  the  ordinate  to  the  axial  distance  of  the  beam  from  the  left  boundary 
normalized  by  the  layer  of  the  beam. 
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PS*] 


Figure  4.2.  Deviation  of  the  Natural  Frequencies  of  a  Concentric  Circular  Conic  Beam  with  a 
Concentric  Circular  Cylindric  Bore  From  those  of  an  Equivalent  Circular  Tube 
as  Function  of  the  Cone  Radii 


max 


1.0 


X/L 

Figure  4.3.  Deflection  of  a  Concentric  Circular  Conic  Beam  with  a  Concentric  Circular  Cylindric 
Bore  Versus  Axial  Distance  for  the  First  Five  Principal  Modes  of  Lateral  Vibration 
for  R  =  .95123 
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